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Abstract. Spiral galaxies that are deficient in neutral Hydrogen are observed on the outskirts of the Virgo cluster. 
If their orbits have crossed the inner parts of the cluster, their interstellar gas may have been lost through ram 
pressure stripping by the hot X-ray emitting gas of the cluster. We estimate the maximum radius out to which 
galaxies can bounce out of a virialized system using analytical arguments and cosmological iV-body simulations. In 
particular, we derive an expression for the turnaround radius in a flat cosmology with a cosmological constant that 
is simpler than previously derived expressions. We find that the maximum radius reached by infalling galaxies as 
they bounce out of their cluster is roughly between 1 and 2.5 virial radii. Comparing to the virial radius of the Virgo 
cluster, which we estimate from X-ray observations, these Hi-deficient galaxies appear to lie significantly further 
away from the cluster center. Therefore, if their distances to the cluster core are correct, the H I-deficient spiral 
galaxies found outside of the Virgo cluster cannot have lost their gas by ram pressure from the hot intracluster 
gas. 
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1. Introduction 

Radio observations at 21cm have revealed that spiral 
gala xies within clusters are deficient in neutra l Hydrogen 
(e.g. IChamaraux. Balkowski. fc Gerard! Il9 80h and their 
H I-deficiency, normalized to their optical diameter and 
morpholog ical type, is largest for the spirals near the clus- 
ter center ilHaynes fc Giovanellilll98rilCayatte et alll990l: 
ISolanes et al.ll200ljh IChamaraux et al.1 l)l98d(l suggested 
that the H I-deficiency of cluster spirals was caused by 
the ram pressure stripping of their interstellar Hydrogen 
by the hot intracluster gas that emits in X-rays. Galaxies 
falling face-on into a cluster experience a ram pressure 
that scales as p c iv 2 l|Gunn fc Gottlll972|h where p c \ is the 
cluster gas density and v is the relative velocity of the spi- 
ral galaxy in its cluster. Therefore, ram pressure stripping 
requires the large infall velocities present in rich clusters. 
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ISolanes et alJ l|2002() recently discovered deficient spi- 
rals in the periphery of the Virgo cluster, with several 
ones typically over 5 Mp c in front or b e hind the cluster 
core. In an ensuing study, I Sanchis et al.l l(2002f) could not 
discard the possibility that some of these galaxies could 
have passed through the cluster core and in the process 
had their interstellar gas swept out by the ram pressure 
caused by the intracluster hot diffuse gas. 

The idea of galaxies beyond the virial radius having 
passed through the m ain body of a cluster in the past 
has been addressed by iBalogh. Navarro, fc Morris! l|2000h 
in the context of the discovery of reduced star forma- 
tion rate s on the outskirts o f clusters in comparison with 
the field llBalogh et al.lll997lh Using cosmological simula- 
tions, IBalogh et al.l analyzed 6 clusters within a sphere of 
2 times their final virial radius and found that 54 ± 20% 
of the particles between r2oo and 2 r2oo (where ^200 is 
the radius where the mean density of the cluster is 200 
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times the critical density) have actually been inside the 
virial radius of the mai n cluster prog enitor at some ear- 
lier time. Unfortunately, iBaloeh et al.l do not provide any 
precision on the maximum distances that such particles 
that were once within a cluster progenitor can move out 
to. Furthermore, one needs to check if 2 r2oo represents a 
sufficient distance for particles bouncing out of virialised 
structures to explain the H I-deficient galaxies on the out- 
skirts of the Virgo cluster. 

In this paper, we ask whether the H I-deficient galax- 
ies on the outskirts of the Virgo cluster have previously 
passed through the core of the cluster, using both ana- 
lytical arguments and the output of cosmological TV-body 
simulations. In Sect. we describe the TV-body simula- 
tions analyzed in this paper. Next, in Sect.|3 we study the 
structure in radial phase space of dark matter halos of the 
simulations. In Sect. 21 we compute the maximum rebound 
radius, both analytically, making use of the turnaround ra- 
dius of cosmological structures, which we compute in an 
appendix, and by studying the structure of our simulated 
halos in radial phase space, as well as analyzing the or- 
bita l evolution of partic l es in t he cosmological simulations 
of lFukushiee fc Makinol l)200lft . In Sect. GO we estimate the 
virial radius and other virial parameters of the Virgo clus- 
ter to permit the estimation of the Virgo rebound radius 
in physical units. We disc uss our results in Se ct. El 

In a companion paper ijSanchis et alJl200^) , we discuss 
in more detail the origin of the H I-deficiency in galaxies on 
the outskirts of the Virgo cluster, by analyzing 2D slices 
of the 4D phase space (right ascension, declination, dis- 
tance and radial velocity) and comparing them with the 
cosmological iV-body simulations used here. 

2. TV-body simulations 

The iV-body simulations used here to find A ^-bodv r eplica 
of t he Virgo cluster w ere carried out by iNininl l|l999l 
see lHatton et alJ l2003j) in the context of the GALICS 
ijHatton et all) hybrid iV-body/semi- analytic model of hi- 
erarchical galaxy formation. Here we are basically inter- 
ested in the density and velocity fields directly traced 
by dark matter particles. The iV-body simulation con- 
tains 256 3 particles of mass 8.3 x 10 9 M Q in a box of 150 
Mpc size and it is run with a softening length amount- 
ing to a spatial resolution of 29 kpc. The simulation 
was run for a flat universe with cosmological parameters 
fl = 0.333, VL A = 0.667, H = 66.7 km s^ 1 Mpc" 1 , and 
as — 0.88. Once the simulation is run, halos of dark matter 
are detected with a 'Friends-of- Friends' (FoF) algorithm 

1 Davis et al . 1985). with a variable linking length such 
that the minimum mass of the FoF groups is 1.65x 10 11 M Q 
(20 particles) at any time step. With this method, over 

2 x 10 4 halos are detected at the final timestep, correspond- 
ing to the present-day (z = 0) Universe. The GALICS halo 
finder does not allow halos within halos, so that a clus- 
ter, to which is assigned a massive halo, cannot contain 
smaller halos within it. 



3. Halo structure in radial phase space 

Fig. ^ shows the radial phase space diagrams, i.e. radial 
velocity vs. radial distance, both relative to the halo cen- 
ter, for 4 massive halos. We have studied the final output 
of the simulation at z = 0. The centers of the isolated 
halos shown in Fig. ^ are provided by the simulation out- 
put and correspond to the barycenter of the FoF groups 
of particles. Radial distances are normalized to the virial 
radius rioo, corresponding to the radius where the mean 
density is 100 times the critical density of the Universe. 
We use rioo instead of r 2 oo because the former represents 
better the virial radius in universes with n on zero cosmo- 
logical constant ijKitavama fe Sutclll996|) . Velocities are 
normalized to the respective circular velocities at r^o- 

Most of the massive halos in the simulation have a 
phase space diagram similar to those for the 4 halos shown 
in Fig. Q once scaled to the virial radius and circular ve- 
locity at the virial radius. In particular, all plots show a 
virialized region for radii smaller than the rioo (although 
this virialization is not perfect, for example see the ex- 
cess of positive velocity particles in the lower left plot, 
presumably caused by a large group that is bouncing out 
of the cluster), and an infalling region with velocity in- 
creasing with radius and asymptotically reaching the (lin- 
ear) Hubble flow. One clearly notices groups or small clus- 
ters of particles in the outer regions (e.g. at r = 3 rioo m 
the upper right plot), which display Fingers of God pat- 
terns in phase space. The material within 1 or 2 rioo that 
bounces out of the cluster should form a pattern symmet- 
ric to the infalling pattern relative to the zero velocity 
line. It is smeared out by numerical two-body relaxation 
(S. Colombi informed us that the rebounding region is 
seen more sharply when simulations are run with increased 
potential softening lengths that reduce the numerical re- 
laxation). The global aspect of these phase space plots 
is similar to those shown by iFukushige fe Makinol (|200lL 
Fig. 21), at various epochs of their cosmological simula- 
tions (run with a standard — tt m = 1, J7a = — CDM 
cosmology). 

4. Maximum rebound radius 

An inspection of Fig. ^ shows that dark matter par- 
ticles beyond the virialized core and outside the in- 
falling/expanding zone of phase space can reach 2 or 2.5 
times the virial radius, but not any further, and moreover 
come in groups of particles which appear to be tidally 
shredded in phase space. In other words, particles that 
cross through the core of a cluster cannot bounce out be- 
yond 2.5 virial radii. 

Now galaxies are not just particles, but arise within 
particle condensations known as dark matter halos, which 
should arise as vertically-elongated (Fingers of God) par- 
ticle condensations in phase space, and therefore ought 
to avoid the fairly sparse regions of phase space where 
the outermost outgoing particles are seen. The open cir- 
cles in Fig. n indicate the halos without galaxies within 
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Fig. 1. 3D radial phase space plots (radial velocity versus radial distance, both relative to the halo center and normal- 
ized to the virial circular velocity and radius, respectively) of dark matter particles in a ACDM cosmological A^-body 
simulation for the four most massive halos at redshift z = 0. The virial radius, mass and circular velocity at virial 
radius are listed on the right corner of each plot. The open and small closed circles are the identified halos in the 
cosmological simulation respectively without and with galaxies within them. 



them in GALICS. The absence of galaxies within halos is 
a feature of GALICS for halos that cross a larger one (its 
galaxies become part of the larger halo). Empty halos can 
also occur in GALICS for isolated halos in which galax- 
ies have not yet had time to form. The empty halos out- 
side the infalling/expanding region do not extend beyond 
l-7rioo (upper left and lower right panels, with the former 
possibly a member of the group at 1.9rioo). In contrast, 



the normal halos (filled circles in Fig. ^) outside the in- 
falling/expanding region do not extend as far from the 
main halo. The positions of the empty circles in Fig. ^ 
therefore suggests that halos crossing the main halo do 
not bounce out further than 1.7rioo- 

This maximum rebound radius is consistent with a 
close inspection of the rig ht panel of Figure 20 from 
iFukushige fc Makinol l)200lh . which shows that the largest 
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rebound radii in one of their cosmological simulations, 
i.e. the largest radii of a particle that has experienced 
at least one pericenter, is 2 Mpc, occurring at the present 
epoch, for a cluster whose present-day virial radius (r2oo) 
is at 1.7 Mpc (see their Table 2). Hence for that partic- 
ular shell, the rebound radius is only 1.2 times r2oo and 
an even smaller factor times r^g (the canonical radius 
for the cosmology used). Given the f2 m = 1 cosmology 
used, for which the spherical infall model yields scale-free 
growth, the rebound radius should be proportional to the 
turnaround radius, which itself should be proportio nal to 
the v irial radius, with a time growth of r ~ i 8 / 9 l|Gottl 
Il975l) . We checked that the other rebound radii occurring 
earlier were even smaller than the scaled expectation of a 
2 Mpc radius today. 

One can confirm this result through simple analytical 
arguments. First, if one identifies the virial radius to the 
radius where infalling shells meet the rebounding shells, 
this will be very close to the rebound radius itself, defined 
as the radius where a shell reaches its second apocenter 
(see Fig. 1 of lMamorJ ri992L so that the rebound radius 
will be very close to rioo- 

Moreover, one can estimate the rebound radius for 
a flat cosmology in the following manner. To begin, as- 
sume that the rebound radius is r times smaller than the 
turnaround radius (the first apocenter of the shell), and 
occurs at a time equal to t times the epoch of turnaround. 
The mass within a shell that has reached its second apoc- 
enter will be close to but greater than the mass M within 
the same shell at turnaround, since some additional mat- 
ter will be infalling for the first time. We write this as 



where the virial radius 



is defined such that the mean 



Preb(^o) Z r 3 p ta (Tta) , 



(1) 



where the rebound density is for the present epoch, while 
the turnaround density is for the epoch of a shell's first 
apocenter given that its second apocenter is today. Given 
that the mass within the shell at turnaround is 

M = y(l + <5i)pirf 
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where we used equations IjA.lOf) and (|A.17|I . We then ob- 
tain a mean density within the turnaround radius that 
satisfies 
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2/ta 
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Equations and lead to 
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where p Cj o is the present day critical density. Inverting 
mean density into mass, assuming din p/dhir = —a, leads 
to 



rvir * I A 



-1/a / ~ N -3/q 
2/ta 



(5) 



density within it is A t imes the critical density. For 
iNavarro. Frenk. & White] l|l995l hereafter NFW) density 
profiles, with c oncen tration parameters c ~ 6, as found by 
iBullor.k et all ll200lh in recent ACDM cosmological simu- 
lations, using A ~ 100, we obtain a slope a = 2.4 at the 
virial radius. 

There are two ways one can reach a large rebound 
radius today: either with a large ratio of rebound to 
turnaround radius for a given shell (i.e. r as small as 
possible, but it cannot be smaller than unity), or with 
a turnaround as late as possible (i.e. t as small as possi- 
ble), which implies a lower turnaround density, hence a 
lower rebound density (for given r). 

In the first case, the rebound radius could be equal 
to the initial turnaround radius, i.e. r — 1. One then ex- 
pects that a shell that is presently at its second apocenter, 
should have reached turnaround at epoch to/3, i.e. t — 3. 
We make use of the reasonably simple expression for the 
turnaround radius versus time obtained in the Appendix. 
Solving equation (|A.21|) for So , using equations (|A. 1 R|> and 
(|A"T9)l for y ta and (O = 0.3, A = 0.7), yields S = 2.69 
and 2/t a = 0.229. Equation (J5J and a = 2.4 then lead to 
r rcb <; 1.8 r ion- 

Second, in the spirit of the shell rapidly virializing 
through violent relaxation at full collapse, one can as- 
sume that the shell reaches its second apocenter right af- 
ter its first pericenter (full collapse), i.e. t ^> 2, with a 
radius near half its turnaround radius (r = 2). This yields 
(not surprisingly) r re b = 1-03 rioo- The very best case re- 
quires r = 1 and t = 2, and leads to r rc b = 2.46 rioo- 
Table ^ summarizes these estimates of the rebound ra- 
dius, where 5q and j/ta come from equations (IA.18(I . (|A.19(1 
and (|A.2ip and r ro b/Vioo from equation JSJ). The last four 
lines are for the four particle orbits from the simulations of 
iFukushige fc Makinol l(200ll) that rebound within the last 
7 Gyr, and for these 5q is computed so that the epoch of 
the second apocenter is the present time, hence the epoch 
of turnaround is to / (t le h / tt a ) ■ 

Table 1. Rebound radius in different scenarios 
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1.0 


3.00 


2.69 


0.229 
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0.258 


0.87 


2.0 


2.00 


2.15 


0.296 


1.03 


1.0 


2.00 


2.15 


0.296 


2.46 


1.5 


2.50 


2.42 


0.258 


1.25 


2.3 


2.50 


2.42 


0.258 


0.73 


2.4 


2.17 


2.24 


0.282 


0.77 


2.3 


3.58 


2.99 


0.205 


0.55 



Although, the uncertain effects of relaxation as the 
shell crosses through the virialized region lead to uncer- 
tain values of the ratios of rebound to turnaround ra- 
dius (r) and time (i) , in the most favorable spherical in- 
fall model, the rebound radius cannot be greater than 2.5 
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times the virial radiu s. The cosmological simulations of 
iFukushiee k, Makinol suggest rebound radii of order of 
rioo, but recall that they are for a different cosmology 
(fl = 1,A = 0). The use of real simulations has the 
added advantage of incorporating the effects of two-body 
encounters that can push material beyond the theoretical 
rebound radius. However, it is not always easy to distin- 
guish in a given snapshot of phase space (e.g. Fig. 0) the 
material that is bouncing out of a structure with mate- 
rial that is infalling for the first time (except that the 
former particles are in halos without galaxies), and the 
analysis of particle h istories, done for the simulations of 
IFukushiee fc Makinol is beyond the scope of this paper 
for the GALICS simulations. 

5. Virial radius, mass and velocity of the Virgo 
cluster 

For the application of results of Sect. 0]to the Virgo clus- 
ter, we require an estimate of the virial radius of the clus- 
ter. We estimate r W o through the large-scale X-ray ob- 
servations of the Virgo cluster obtained with the ROSAT 
All-Sky Survey bv lSchindler. Binggeli. fc Bohringerll|l999l) . 
The peak of the X-ray emission in Virgo coincides with 
the position of the giant elliptical M87, and we use 
the integrated mass profile around M87 obtained by 
ISchindler et alJ to derive the virial radius. With their 
isothermal approximation for the total mass profile (their 
Fig. 11a), M(r)/r is independent of radius for r ~j>> r c and 
also of the assumed distance to Virgo, yielding 



GMW =3p kT 



(G) 



where (3 = 0.47 is the shape parameter of the X-ray gas 
density profile, kT = 2.7 keV and iim p is the mean particle 
mass, generally assumed for a hot plasma to be roughly 
0.6 uip-t where m p is the proton mass. The virial radius, 
rioo, is then obtained through 



3M 



ion 



A 3jff ° 2 



;ttG 
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where A = 100 is the spherical overdensity at the virial 
radius. Solving equations © and (JJJ), one obtains 



T-ioo 



M 



100 



wioo 




= 1.65 /i"^ Mpc , 



/jLTTip 

2.30 x 10 14 h~l M 







G M 



100 



7*100 



l3(3kT 



780kms~ 



(8) 

(9) 
(10) 



where i>ioo is the circular velocity at rioo and h 2 /3 = 
#o/(66.7 km s^Mpc" 1 ). 

The Virgo cluster is believed to have a complex 
structure, as it may be composed of several subclusters 



l|Binggeli. Tammann. fc Sandagell987)) around the ellipti- 
cal galaxies M87 (NGC 4486), M86 (NGC 4406), and M49 
(NGC 4472). The virial radius of 1.65 h~LMpc is such 
that the important substructure surrounding M86 is well 
within it, at a distance of 0.23 (projected) and 0.44 (3D) 
times rioo (converted to H = 70kms~ 1 Mpc -1 ) of M87. 
Similarly, the important substructure around M49 lies just 
within the virial radius, at a distance of 0.82 (projected) 
and 0.88 (3D) times the virial radius of M87. Among the 
other Messier galaxies in Virgo, which all have fairly se- 
cure distances, M60 is well within the virial radius, M89 
lies just at the virial radius, while M59 is outside. 

Analyzing the X-r ay emission of the intra-cluster gas, 
ISchindler et all l)l999h found that the M49 subcluster is 
2.4 times less massive than the M87 subcluster. The ratio 
of masses b etween the M87 and t he M86 subclusters is 
even larger l|Bohringer et al.lll994h . It therefore appears 
that the structure surrounding M87 is by far the most 
massive component in the Virgo cluster, and it is not a 
bad approximation to choose a single halo to represent 
the cluster, so we can apply the results of Sect. 0] 



6. Discussion 

The analysis of Sect. 0] indicates that the maximum re- 
bound radius is between 1 and 2.5 times the virial radius. 
Given the virial radius of 1.65 Mpc for the Virgo clus- 
ter, which we derived in Sect. and the maximum re- 
bound radius derived in Sect. 01 galaxies passing through 
the Virgo cluster core in the past cannot lie further than 
(1 — 2.5) x 1.65 = 1.7 — 4.1 Mpc from the cluster center. 

iBalogh et all (2000) found that a very significant frac- 
tion of particles at a distance between 1 and 2 r2oo from the 
centres of cosmologically simulated clusters have passed 
through the main body (r < r2oo) of a cluster progeni- 
tor at some earlier epoch. Note that although the analy- 
sis of IBalogh et alJ was performed in terms of r2oo, they 
would have very probably gotten similar results had they 
scaled their clusters and progenitors with rioo instead. As 
mentioned in Sect. ^ it is not clear from their analysis 
if particles can escape beyond 2 r2oo or even beyond, say, 
only 1.5r 2 oo- Also, it is easier to displace out to large dis- 
tances particles rather than large groups of particles rep- 
re senting a gala xy (or subhalo). In any event, the result 
of IBalogh et alJ is consistent wit h our analysis of Sect ,0] 

An examination of Fig. 2 of ISolanes et al.l l)2002h in- 
dicates 3(7 Hi-deficient galaxies lying between 9 and 30 
Mpc from the Local Group, and in particular galaxies at 
10 and 28 Mpc from the Local Group, whose distance error 
bars do not reach the wide range of distan ces to the Virgo 
cluste r found in th e literature (14 Mpc bv lGiardullo et alJ 
ll99Sl to 21 Mpc bv lEkholm et alJl200(t . Therefore, it ap- 
pears very difficult to explain such H I-deficient galaxies 
over 5 Mpc in front or behind the cluster center as hav- 
ing crossed through the center of the cluster and bounced 
out if their distance estimates are accurate. This would 
suggest that the H I-deficient galaxies in the outskirts of 
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the Virgo have not had their interstellar gas ram pressure 
stripped by the intracluster diffuse hot gas. 

An alternative explanation to the presence on the out- 
skirts of clusters of Hi-deficient spirals, as well as to the 
decreased star formation rates and redder colours of galax- 
ies in these regions, relative to field galaxies, is that the 
three effects of H I removal, decreased star formation and 
redder colours, all intimately linked, may be caused by a 
significant enhancement of massive groups of galaxies at 
the outskirts of clusters, as exp ected from t he statistics 
of the primordial density field feaiserlll984 | )_ applied to 
small groups versus rich clusters l|Mamonlll995|) . If this 
is the case, we would then expect a correlation between 
H I-deficiency and X-ray emission from the intragroup gas. 
However, while tidal effects, which to first order depend 
on m ean density, regardless of orbit eccentricity ijMamonl 
l200(t . are similar between less massive groups and more 
massive clusters, ram pressure stripping effects, which also 
depend on the squared velocity dispersion of the environ- 
men t will be muc h reduced in groups relative to clusters 
fe.g.lAbadi et aklluffll. 

In the companion paper l|Sanchis et alJl2003|) . we con- 
sider different explanations to the origin of the Hl- 
deficiency of these outlying galaxies: 1) incorrect dis- 
tances, so these objects would in fact lie close enough to 
the cluster core to be within the rebound radius and their 
gas could have been removed by ram pressure stripping, 
2) incorrect estimation of the H I-deficiencies and 3) tidal 
perturbations (stripping or heating) by nearby compan- 
ions or within groups. 
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Appendix A: Radius, time and density of 
turnaround in a flat ACDM cosmology 

In this appendix, we compute the parameters of shells at 
turnaround, i.e. reaching their first apocenter, in a ACDM 
Universe without quintessence (wq = — 1). 

In a non-quintessential Universe with a cosmological 
constant, the equation of motion of a shell of matter is 



d 2 R GM(R, t) A 



dt 2 



R 2 



3 



R 



(A.l) 



where the first term on the right hand side is the gravita- 
tional force and the second term is the repulsive effect of 
the cosmological constant A. A given shell of matter that 
first expands with the Universe and then turns around 



and collapses will begin to cross shells that have already 
settled in a structure only slightly before its own collapse. 
Therefore, during the expansion phase, there is no shell 
crossing and hence M(R,t) = est. 

The equation of motion can then be easily integrated 
to yield the energy equation 



E = I f z - 



GM 1 



1 



- -Ar z = - 
6 



GM 



l^r 2 , (A.2) 



where E is the energy per unit mass of the shell, and where 
the i subscript refers to a very early time U , corresponding 
to redshift 3> 1 (e.g. Zi — 1000). It is convenient to use 
the quantities 



8nGp 

Yh 2 



and 



A 



A 
3H 2 



(A.3) 



which represent the dimensionless mass density and dark 
energy of the Universe at any epoch, and we will use sub- 
scripts '0' to denote the present epoch (z = 0). In the case 
of a flat Universe (f2rj + Ao = 1), as confirmed w ith the 
WMAP CMB experiment bv ISoereel et alJ l|2003|) . mass 
and energy conservation applied to t he Universe lead to 
(e.g. iRichstone. Loeh. fc Tumerlll99l 



Hf =H 2 [fi (l + 2 4 ) 3 + A ] 



(A.4) 



At time h, the mass enclosed within radius r,- is 



M = ^ (1 + Si) Sk Pc:i rf = \ (1 + Si) Qi H 2 1± , (A.5) 

where Si is the relative overdensity within radius Vi at 
time ti (normalized to the density of the Universe at that 
epoch), and where we made use of the critical density of 
the Universe at L: 



3HI 
8ttG 



(A.6) 



With a dimensionless radial growth factor u(r,t) 
r/ri, the first energy equation can be expressed as 



. TT 2E (i + <yfii , _ 

u=Hi \m + — u — +XiU ' 



(A.7) 



where we made use of equations (|A.2ll . (|A.3|) . (|A.4|I . and 

(JA.5(1 . The energy E of the shell is obtained by expressing 
equation l|A.7|) for the epoch ti, yielding 

7^ = SJ-Kl + W + Ai] • (A.8) 

i i i 

The second equality of equation (|A.3|I yields 

2 /r-> \ -l 



A; 



0((1 



< 1 



(A.9) 



Mass conservation amounts to a mean density of the 
Universe that varies as 



Pi = p(zi) = po (1 + ZiY 



(A.10) 
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where po is the present-day mass density of the Universe. The time of turnaround is then obtained by integrating 
This then yields equation l|A.14|) and writing u = (1 + Zi) y, yielding 



An Q l 1 The starting point of the integration is taken at instead 

— 1 — ^' (A-ll) Q £ _|_ which corresponds to adding the negligible 



time between the t = and U. With equation (|A.4() . eaua- 
Inserting equations (|A.8|) . (| A. 9|> and (| A. 1 1|> into equa- ^ion ||A.20|l becomes 
tion yields 

^oT ta = f 5 5 =* , • (A.21) 



U = H *\TFT^ 1- A, - (1 + dj) , (A. 12) 

V H i u For fi = 1, A = 0, we get y ta = 3/(5 S ), H T ta = 

, j j , » s \ \ ~\i ■ i j j (7r/2) [3/(5(5o)1 3//2 , so that a shell that collapses today, 

w here we discarded ter ms m o(A 7 ;). Now, as pointed out v ' ' 1 ' v U/J ' , J ' 

bvlChodorowskiim.lBartehnann. Ehlers. fc Schneider! hence tu ™ s ™ d at ^°^ a = ^ *°/ 2 = \ /3 ! re " 

nqq^ anHlPaHmflnflbiiflnliliqq^ nnP qbnnlH not aTOlm p quires a lmearl y extrapolated density contrast of S = 

^liiijct) aiui i.r aoiiiaiiatjiiani iji^hMii . one snouici not assume ("3/5) (3 /2)^/ 3 1 686 t d 

that the initial flow is a pure Hubble flow (u, = if,), \ I )\ I ) ~ > P 

but incor porate the peculi ar motion acqui red before Zj. 

Usi ng the IZeldovichl jl97Cl) approximation. IChodorowskil References 

and| Bartelmanne laL|show that Abadi; M q Moor6; g & BowCT; R q 1999; M NRAS, 

i. = Wl - ^ (A 13) 308 ' 947 

1 ^ 3 J ' V } Balogh, M. L., Morris, S. L., Yec, H. K. C, Carlberg, 

__, . . , R. G., & Ellingson, E. 1997, ApJ, 488, L75 

Inserting equation 1A.13I) into equation (IA.12I) and dis- „ , , •., T , T t e i\/r • c t onnn a t 

° . 7 ' \ ^ , -, r— ■ — ' . . Balogh, M. L., Navarro, J. F., & Morris, b. L. 2000, ApJ, 

carding high order terms m (1 + Zi) and \i — writing 

Si = 0(1/(1 + z^) = 6 /(l + z^ — yields Bartelmann, M., Ehlers, J., & Schneider, P. 1993, A&A, 

j—^ 1 280, 351 

u = Hi J-~ Si + - + A, u 2 . (A. 14) Binggeli, B., Tammann, G. A., & Sandage, A. 1987, AJ, 

94, 251 

The turnaround radius is obtained by solving for u = Bohringer, H., Briel, U. G., Schwarz, R. A, et al. 1994, 

0, i.e. solving Nature, 368, 828 

3 5 Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, 

*i«ta- g*i«ta + l = (A.15) MNRAS, 321, 559 

e m 9 tj a ira ( j i-i , ,. Cayatte, V., van Gorkom, J. H., Balkowski, C, & Kotanyi, 

for Mt„. ILokas & Hoffman! (12001 ) found the same solution ^ 1990 AJ 100 604 
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We can go further than| Lokas fc Hoffma,n| by expressing cho d orowski , M. 1988, Master's thesis, University of 

and solving a cubic equation in terms of present-day quan- Warsaw 
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solvmg R. E. 1998, ApJ, 492, 62 
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r =r i (l + z l ) (A.17) Fukushige, T. & Makino, J. 2001, ApJ, 557, 533 
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is the comoving radius of the initial perturbatio n. The Qunn, J. E. & Gott, J. R. 1972, ApJ, 176, 1 

smallest real positive solution of the cubic equation HAT6J) Hatton, S., Devriendt, J., Ninin, S., et al. 2003, MNRAS, 

is 343, 75 

2 y/5 I S Q n Q f(/) + Tr\ Haynes, M. P. & Giovanelli, R. 1986, ApJ, 306, 466 
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(l-fio)M, 



t 500 
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